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Q Abstract 

We review various inequalities for Mills' ratio (1 — $)/</>, where 4> an d denote 
the standard Gaussian density and distribution function, respectively. Elementary 
considerations involving finite continued fractions lead to a general approximation 
i ' scheme which implies and refines several known bounds. 

1 Introduction 

a 

Explicit formulae for the distribution function $ of the standard Gaussian distribu- 
te i 

>• tion are unknown, apart from various expansions, e.g. the series expansion 

m 

1 1 ~ (_1\k T 2k+l 

cm' k ~° 

for real x, or the continued fractions expansion 

(f)(x) 



5-H 



1 - $(x) 



X x + 



& x+ 3" 

x H 



2 



for a; > 0. Here </> = $' denotes the standard Gaussian density. We refer to 
Abramowitz and Stegun (1972, Chapter 7) for these and numerous further results 
about the function $. The expansion ([I]) indicates that there should be good ap- 
proximations or bounds on 1 — of the form 

!-$(*) » ^ *>0, 
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with h : (0, oo) — > (0, oo) being a relatively simple function. For instance, let 
ho(x) := x, hi(x) := x + 1/x and, for k > 2, 

h k (x) : = xH . 

x H 



x H — 

x 



Then it is known that for integers k > and arbitrary x > 0, 

, , j < 1 <K X ) -r ; • I even ' 

!>} sM s odd. 



In particular, 



< 1 - $(x) < 



x + 1/x ' X ' 

which was first established by Gordon (1941). These special functions and their 
properties have been investigated by numerous authors; we only refer to Shenton 
(1954), Pinelis (2002), Baricz (2008) and the references therein. 

While the previous bounds are only useful for x bounded away from zero, some 
authors provided inequalities on the whole interval [0, oo) or even larger sets. Indeed, 
Komatu (1955) showed that 

26(x) , s 26(x) 

(3) . yV 7 < l-$(x) < . yV ; forx>0; 

V4 + x 2 + x V2 + x 2 + x 

see also Ito and McKean (1974). In fact, the lower bound is due to Birnbaum (1942) 
who formulated it equivalently as 



a/4 + x 2 — x 



(x) < 1 - $(x) for x > 0. 
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Pollak (1956) refined Komatu's upper bound as follows: 

(4) l-$(x) < ; 2 ^ forx>0. 

v/8/tt + x 2 + x 

An alternative upper bound, due to Sampford (1953) and rediscovered by Szarek 
and Werner (1999), reads 

(5) l-$(x) < , forx>0. 

V8Tx^ + 3x 

Shenton (1954) and Kouba (2006) generalized Komatu's bound ^ and Sampford's 
bound ^ substantially. Here is a reformulation of their bounds with continued 
fractions: For integers k > and j = 1, 2 define h^j : [0, oo) —> (0, oo) via 

h Qtj {x) := y/j/2 + (x/2) 2 + x/2, h ltj (x) := x+ ' 



a /1+j72 + (x/2) 2 + x/2' 



and, if k > 2, 

h kJ (x) = x + 



x H 

x + 



^k + j/2 + (x/2) 2 + x/2 
Then for all x > and integers m > 0, 

(6) < 1 _ < 



h 2m ,2( x ) h 2m ,i(x) 

/„\ 0(^) -, ^/ N 

(7) / x < 1 - < 



^2m+l,l(^) ^2m+l,2(^) 

In the present manuscript we present all these bounds in a common framework 
and propose refinements. In Section[2]we consider the derivative of (j)/h — (1 — $) for 
a smooth function /i : (0, oo) — > (0, oo) with lim^oo = oo, and these elementary 
considerations yield the bounds ([3]), @, ([5]) and a new lower bound. In Section [3] 
we consider approximations <f)/h k of 1 — $, where 

1 



h k (x) 



x + 



k-l 

x H 



A; 

x + 



<7*W 

with smooth functions g k : (0, oo) — >■ (0, oo). It turns out that under general condi- 
tions on g k , this yields upper and lower bounds <fi/h k of 1 — $. 

In Section |4j we consider functions g k {x) of the form 

g k {x) = y/c k + {x/2) 2 + x/2 

with special constants c k E [k + 1/2, k + 1], and this improves the bounds ^ and 
(J7|). In Section [5] we consider 

9k{x) = y/Ck~+X k X 

with an additional constant X k > 0, which leads to purely rational functions h k . 
The resulting bounds are compared to those in Section [4} These rational functions 
h k may be improved substantially ba considering 



gk(x) = y/ck~exp(-5 k x) + x 

for some 5 k > as explained in Section [6j 

Finally, in Section [7] we describe briefly a recurrence scheme to represent the 
functions h k as ratios p k /q k rather than continued fractions. 



3 



2 First Steps 



Let h : (0, oo) — > (0, oo) be a differentiable function with linx^oo h(x) = oo. Then 
the approximation error 

A <*> = fg} - (i - 

satisfies 

(8) lim A(x) = 

X— ¥oo 

and 

(9) A'(x) = ^(h(x) 2 -xh(x)-h'(x)). 

Consider first h(x) = x. Then h(x) 2 — xh(x) — h'(x) = —1, so A' < on (0, oo). 
Together with (Jsj) we obtain that A > 0, which is inequality ^ for k = 0. However, 
in view of the bounds ^ and Q of Komatu and Pollak, we try 



a/4c + x 2 + 



.r 



:= (x/2) 2 + x/2 

for some c > to be specified later. Then A is well defined and continuous on 
[0, oo), and one verifies easily that 

h(xf - xh{x) - h'{x) = c- 1/2 -t(x,c)/2 

where 

t(x, c) 



V4:C + x 2 ' 

Note that t(-, c) : [0, oo) — > [0, 1) is bijective and increasing. In case of c = 1, 

sign(A'(x)) = sign(l — t(x, c)) = +1. 

Thus A < on [0,oo), and we rediscover Komatu's lower bound in Setting 
c = 1/2, we realize that 

sign(A'(ic)) = sign(— t(x, c)) = —1, 

so A > on [0, oo), which implies the upper bound in 
Note that 

A(0) = if, and only if, h(0) = y/2/n. 



4 



The latter condition is satisfied if c = 2/tt, which corresponds to Pollak's function 



v/8/tt + x 2 + x 
h{x) = . 

Indeed, 

sign(A'(x)) = sign(2c — 1 — t(x, c)) = sign(^ t(x,2/n)j. 

Thus A' > on (0, x a ) and A' < on (x , oo), where x solves the equation 
t(x ,2/ir) = (4 — 7r)/7r G (0,1). This shows that A > on (0,oo), and we ob- 
tain Pollak's upper bound Q. 

Now we go one step further: Since 

\jAc + x 2 + x \/Ac + x 2 — x 2c 
x -\ = x + 



2 2 y/*c + x 2 + x' 

we consider functions h of the form 

1 

h[x) = x + 



with g : (0, oo) — > (0, oo) differentiable. One can easily verify that 

h[x) — xh[x) — h \x) = — — , 

g(x) 2 

so 

sign(A'(x)) = sign(/i(x) 2 — xh(x) — h'{x)) 

= -sign^x) 2 - xg(x) - 1 - g'(x)). 

In case of g(x) = x, we obtain sign(A'(x)) = +1, so A > on (0, oo), i.e. (|2j) 
holds for k — 1. Again we can refine this considerably by considering 



g{x) = + (x/2) 2 + x/2 

for some c > 0. This leads to 

, . . 2 V4c + x 2 - x 

h[x) = x H . = x H 

\/4c + x 2 + x 2c 

(2c- l)x + V4c + x 2 

2^ ' 

and 

sign(A'(x)) = — sign(c — 1 — g'(x)) = — sign(2c — 3 — t(x, c)) . 
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For c = 2, the latter sign equals —1 for all x > 0, so A > on [0, oo), and we 
obtain the upper bound ^ of Sampford. On the other hand, if c = 7r/2, then 
h(0) = 1/0(0) = so A(0) = 0, and 

sign(A'(x)) = — sign(7r — 3 — t(x, vr/2)). 

Thus A' < on (0, x a ) and A' > on (x a , 00), where t(x ,7t/2) = n - 3 G (0,1), 
and we obtain a new lower bound: 

(10) 1 - $(x) > for x > 0. 

(-7T — l)x + V27T + X 2 

This corresponds to 

2 

h(x) = x H — . 

V2vr + x 2 + x 

which is strictly smaller than the function 

2 

ri(x) = x + 



a/4 + x 2 + x 
corresponding to Komatu's lower bound in (|3]). 

Figure [T] illustrates the bounds we have seen so far. Precisely, it shows the 
approximation error A = (fi/h — (1 — $) for Komatu's lower bound ([3]), the new 
lower bound ( 10 ) as well as the upper bounds Q of Pollak and (|5| of Sampford. 



3 Continued Fractions and General Bounds 

Recall that we started with an arbitrary function h(x), then turned our attention 
to h(x) = x + l/g(x) with smooth g, and the special functions g we used may also 
be written as g(x) = x + l/g(x) with another smooth function g > 0. After playing 
around with the resulting approximation error A = <p/h — (1 — $) and sign(A'), the 
following scheme seems to be promising: 

h( x ) = 9o(%) = — = x^ 

9i(x) x + —— 

for some differentiable function g 2 : (0, 00) — > (0, 00). Indeed, elementary calcula- 
tions reveal that 

sign (A' (x)) = sign (5-0 (x) 2 - xh (x) - - g' (x)) 
= - signal (x) 2 -xg x (x) - 1 - g[(x)) 
= sign(g 2 (x) 2 - xg 2 (x) - 2 - g' 2 (x)), 

and this suggests a more general result which will be proved via induction: 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Figure 1: Approximation errors A after first steps. 



Lemma 1. Let go, gi,g 2 ,... be differentiable functions from (0, oo) to (0, oo), and 
define 

h (x) := g (x), hi_(x) := xH ]—, h 2 {x) : = xH ^ — 

9i (x) 



X ~\~ 

92{X) 

and, for integers k > 2, 

h k (x) := x H . 



x + 



k-1 

x H 



k 

x + 



9k{x) 

Then the approximation errors A k := <j)/h k — (f — $) satisfy 

s ign(A' k (x)) = (~ l ) k sign(9k(x) 2 - xg k (x) - k -g' k (x)). 

Proof of Lemma [TJ We know that the claim is correct for k — 0,1. Now suppose 
that it is correct for an arbitrary integer k > 1, and let g k (x) = x + {k + l)/g k+1 (x) 



for a differentiable function g k+ \ : (0, oo) — > (0, oo). Then 

g k (xf - xg k (x) -k- g k (x) 

= 9k(x){g k (x) -x) -k-g' k (x) 

fe + 1 \ fc + 1 (H1 K +I (x) 
x H -— -— — k — 1 H 



fl f *+iW / fl , *+iW ^fc+i^) 2 

(#fc+i(£) 2 ~^ft+i(^) - + !) -fljfc+iW)' 



5fc+i(^) 2 



In particular, sign (A' (x)) equals 

signU fc+ i(x) 2 - xg k+1 (x) - (k + 1) - 



□ 

With Lemma [T] at hand we can derive bounds for 1 — $, similarly as in Section [2] 
Note first that g k (x) := x yields g k (x) 2 — xg k (x) — k — g k (x) = —(k + 1), so the 
corresponding approximation error A satisfies sign(A'(x)) = (— l) fc+1 . This implies 
^ for arbitrary k > 1. But Lemma [I] leads to a refined criterion: 

Lemma 2. In the setting of Lemma [IJ suppose that g k is defined and continuous 
on [0,oo) such that h k (0) = v 2/7T, i.e. Afc(0) = 0. If there exists a point x k > 
suci tiat 

sign(#fc(x) 2 - xg k (x) - k - g' k (x)) = sign(x fc -x), 
then A k > on (0, oo) for even k, and A k < on (0, oo) for odd k. 

The requirement that fofc(0) = a/2/tt for fc = 0, 1, 2, . . . means that c£ := #fc(0) 2 
satisfies 

2 A; 2 

(11) Cq = — and c* k = — — for k = 1, 2, 3, . . . . 

^ c k-i 

In other words, 

(12) c* = ^ ct = | and 4 = (^A^'c*^ for k = 2, 3, 4, . . . . 



This leads to 
(13) 



2-4-6 •■■ (2m) \2 2 
2 "' ~ Vl-3-5 ■■■ (2m -1) J ^ ! 



,1.3-5---(2m + l)x27r 

( 14 ) c 2m+1 



2 • 4 • • • (2m) 7 2 

for integers m > 1. For later purposes it is crucial to have good bounds for these 
constants c* k . Numerical experiments led to the formulation of the following result: 
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Lemma 3. The constants c* k just introduced satisfy 



1 



8(k + l] 



< c* k - k - 1/2 < 



1 



< 



1 



\(k + l/2) 



for all k > 0. 



The proof of this lemma will be postponed to the next section, because there we 
get the initial bound k + l/2<c* k <k + l almost for free. 



Remark 1. Note that for even integers k > 2, 



k/2 



so Lemma [3] yields 



k 



k/2 )2-^ 1 ^^kTl/2 



<k + 1/2 



16(fc + l/2) 2: 



4 Refining Shenton's Bounds 

Starting from Lemma [T] we consider 

9k(x) := \/c k + (x/2) 2 + x/2 > x 
for some constant c k > yet to be specified. Now 

g k {xf - xg k (x) - k - g' k (x) = c k - k - 1/2 - t(x, c k )/2. 

Thus 

sign(A' fc (x)) = (-l) fc sign(2(c fc -A;-l/2)-t(x,c fc )). 

In case of c k — k + 1, the right hand side equals (— l) fe sign(l — t(x,c k )) = (— l) k , 
and for c k = k + 1/2, the right hand side equals (— l) k sign(— t(x, c k )) = (— l) fe+1 
for all x > 0. This leads to Shenton's and Kouba's strikt bounds in ^ and ([7]). It 
entails also that the constants c* k from the previous section satisfy 

c* G (/c + l/2,fc + l). 

For /ifc(0) is a continuous and strictly monotone function of c k G [k + 1/2, k + 1] 
with extremal values being strictly smaller and strictly larger than \/2/V. Setting 
Cfc = c* k yields a function g k satisfying the criterion of Lemma [2] with x k solving the 
equation 

t(x k ,c* k ) = 2(c* -k -1/2) G (0,1), 
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i.e. 

fl5 , = 2y^(c£ - fc - 1/2) = 2^(4-/^-1/2) 

Xk y/lfi ~ K - k - 1/2)2 ^-^(fc + l-^) 
These considerations yield already the first part of our main result. 

Theorem 1. Let 

h (x) := v / 2/n + (x/2) 2 + x/2, h^x) := x- 



'ir/2 + (x/2) 2 + x/2 
and, for integers k > 2, 

1 



h k (x) := x + 



x + 



x 



yfc + (x/2) 2 + X/2 

Then the approximation errors A f. := 0/hfc — (1 — $) satisfy the following inequalities: 

A > A 2 > A 4 > ■ ■ • > and A x < A 3 < A 5 < - - - < 
on (0, oo). Moreover, 

i 

max |A fc (x)| = |A fc (x fc )| < 



x>o 1 " v 32(£; + l/2) 2 



with x^ given by (15). 

Remark 2. The first three bounds 4>/h k from Theorem [I] are given by the following 
functions h k : 

, / \ a/8/7T + X 2 + X 

h (a;) = 



hi(x) = x + 



2 

1 (tt - l)x + V2tt + x 2 



^/tt/2 + (x/2) 2 + x/2 vr 

1 8/tt 



h 2 (x) = x H s = x + 

(8/tt - l)x + v/32/vr + x 2 



x 



^8/tt + (x/2) 2 + x/2 

Note that the upper bound <f)/h 2 for 1 — $ is better than Pollak's upper bound 0/ho 
in Q. It is also better than the upper bound (j5]) of Sampford, because the latter 
equals <f>/h with 

h(x) = xH = = xH f; < h 2 (x) 

K) y/2 + (x/2) 2 + x/2 x + 2 2 ^ 

v/2 + (x/2) 2 + x/2 

for x > 0. 
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Figure [2] shows the approximation errors A k from Theorem [T] for some choices 
of k. On the left hand side one sees A , Ai, . . . , A 5 . Note that the worst functions 
A , Ai correspond to the bounds Q and ( 10 ), respectively, also depicted in Figure[TJ 
On the right hand side one sees A 2 , A 3 , . . . , A 9 . 
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Figure 2: Approximation errors A k for k = 0, 1, . . . , 5 (left) and k = 2, 3, . . . , 9 
(right). 

Before proving Theorem [T] we have to prove Lemma [3] 



Proof of Lemma [3j As we have just shown, k + 1/2 < c k < k + 1. Now we 
compare c* k with d k '■— k + 1/2 + (8(k + 7)) 1 for some fixed 7 > 0. Note that 



(16) 

Suppose we can show that 
(17) 



lim ^ = lim ^ = 1. 

k— >oo k k— >oo k 



< £*±i for all A: > 0. 

dk dk+2 



Then, for arbitrary integers k > 0, the sequence (c k +2i/ 'dk+u)iLo is strictly increasing 
with limit 

1- C l+2l y C *k+2il + 2^) 

hm = lim — — = 1, 

£->oo d k+2 e 6-^00 d k+ 2£/ [k + 2£) 

whence c* k < d k . Analogously, if 
(18) 



c * > £*±i for all A; > 0, 



fc+2 



then c* k > d k for each k > 0. 
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Writing x ~ y if x = ry with r > 0, we may deduce from (12) that 



C k+2 _ ^fc 
dk+2 dk 



(k + 2fd k - (k + l) 2 4 +2 
(k + 2) 2 4 — (k + l) 2 d k+2 
(k + 2) 2 (k + 1/2) — (k + l) 2 (k + 5/2) + 



+ 



(fc + 2) 2 

8(A; + 7) 8(A; + 2 + 7) 



+ 



1 (fe + 2) 2 

2 8(A; + 7) 8(A; + 2 + 7) 



(& + 2) 2 (A; + 2 + 7) - 4(fc + 7) (A; + 2 + 7) - (k + l) 2 (k + 7) 
(3 - 6-f)k + 8 - 57 - 4 7 2 . 



In case of 7 = 1/2, the previous expression equals 9/2, so (17) is satisfied. In case 



of 7 = 1 we get —3k — 1, so (18) holds true. 



It remains to prove the refined upper bound c* k < d k := k + 1/2 + (8c* k ) for 



arbitrary k > 0. Since the latter constants satisfy (16), too, it suffices to verify (17). 



Tedious but elementary manipulations show that this is now equivalent to 

k + 7/4 



c k < k + 1/2 + 



{k + 2) 



and the right hand side is easily shown to be larger than our preliminary bound 
A; + 1/2+ (8(^ + 1/2))^ for c* k . □ 

Proof of Theorem [H Our previous considerations show already that 1—$ < 0//i& 
for even k and 1 — $ > <p/hk for odd k on (0, 00). To verify that (j)/ho > 0//i2 > 
0//i4 > • • • and (f>/hi < 4>/hz < 4>/h§ < • • • , it suffices to show that for any integer 
k > and x > 0, 



^/c*. + (x/2) 2 + x/2 = x + 

is strictly smaller than 
(19) x+ 



^c* k + (x/2)2 + x/2 



fc + 1 



X + 



A; + 2 



^ +2 + (x/2)2 + x/2 



To this end, recall that c£ +2 = c*J p% with p fc := (fc + l)/(k + 2) e [1/2, 1) by (12). 



12 



Consequently, (19) equals 

fe + 1 

x H 



k + 1 

x + 



x + 



v/c* + {p k x/2f + p k x/2 
(k + l)c* k 



c* k x + (k + l) (v 7 ^ + (pfcx/2) 2 - p k x/2) 
c| 

^ + V c l + (P^/ 2 ) 2 + (2c^/(fc + 1) - Pk)x/2 
Hence we have to show that 

y/cl + {p k x/2f + (2c* k / {k + 1) - Pk )x/2 < ^cl + (x/2f + x/2, 

that means, 

(2ct/(k + 1) - p k - l)x/2 < ^cl + {x/2f-^cl + { Pk x/2)^ 

{l-p 2 k ){x/2f 
v /c* + (x/2)2 + ^ci + {p k x/2y 

Dividing both sides by x/2, one realizes that the previous inequality holds for all 
x > if, and only if, 2c* k / {k + 1) — p k — 1 < 0, which is equivalent to 

4 < (<: + 1 , )( \ +3/2) = k + l,2 + 1 



k + 2 ' 2k + 4 

But this inequality is weaker than and thus a consequence of the upper bound in 
Lemma HI 

Let A k := <p/h k — (1 — $), where h k is defined as h k with k + 1/2 in place of 
tip,, so sign(A' fc ) = (— l) fc+1 . Thus < A k < A k with A k strictly decreasing, if k is 
even, while > A^ > A k with A k strictly increasing, if k is odd. In case of even k, 
maXa;>o A k (x) is strictly smaller than 



while in case of odd k, maxj.> |A^(x)| is strictly smaller than 



-A fe (0) = A fe (0)-A fe (0) = ^L(i-^h = ^-J*±l!l\ 



If c 



< ^( 4/^-^-1 



2 v y k + 1 /2 

Since c* k < k + 1/2 + (S{k + 1/2))" 1 = (k + 1/2) (l + (8(k + 1/2) 2 )" 1 ) by Lemma^ 



2\\k + l/2 / 4 VA: + 1/2 / ~ 32(£; + l/2) 2 
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5 Rational Bounds 



It is also possible to obtain approximations (j)/h k of 1 — $ with rational functions 
h k ■ [0,oo) — > (0,oo). To this end, consider h k as in Lemma [l] with the simpler 
function 

9k{x) := ^/c* k + \ k x 
for some Afc G [0, 1] yet to be specified. Now 

g k (x) 2 - xg k (x) — k — g' k (x) 
= ^/cj{2\ k - l)x - A fe (l - X k )x 2 + c* k -k- X k . 

Suppose first that X k — 1. Then the previous display equals \[(? k x — (k + 1 — c* k ), 



so 



. . . , k + l — dt 

sif 



ign(A' fc (x)) = (-l) fc sign(z- 



Hence for even k, <p/h k is a lower bound 1 — $, whereas for odd k it is an upper 
bound. Numerical experiments showed, however, that the bounds in Section [| are 
better. 

More interesting is the choice A^ := c* k — k £ (1/2, 1). Then 
g k {xf - xg k (x) - k - g' k (x) = x(a k - b k x) 



with a k := 2y r c* k {c* k - k - 1/2) and b k := (c* k - k)(k + 1 - c* k ). Thus g k satisfies the 
criterion in Lemma [2] with x k equal to 

i\f^M - k - 1/2) 2Je k {c* k - k - 1/2) 



(20) Xk ' (cj -*)(* + 1-cJ) 1/4- (cj-fc- 1/2)^ 

These considerations yield already the first part of our second main result: 

Theorem 2. Let 

h {x) := y/2jrr + (2/ir)x, ^i(x) := x+ ) 

+ (7T/2 - l)x 

and, for integers k > 2, 

h k (x) := x^^ . 

x H 



a; 



'c* + (c* - k)x 

Then the approximation errors A k := (f)/h k —(l—$) satisfy the following inequalities: 
A > A 2 > A 4 > ■ ■ • > and A\ < A 3 < A 5 < ■ • ■ < 
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on (0, oo). Moreover, the approximation errors A k and A k (as in Theorem satisfy 
the following inequalities: 

max\A k (x)\ = \A k (x k )\ < |A fc (a; fc )| = max I A k (x) I 

x>0 II II I „.^ n I 



x>0 



with x k and x k given by (15) and (20), respectively. But 



Aft 0*0 > |Afc(x)| for x > x k 

( 1 

2xfc < Xfc < min 1, — =7 

V A /ct(l- 



(l-(44)' 2 )^ 



and 



The approximation errors A k are depicted in Figure |3j They look similar to the 
errors A^ in Section^ For a direct comparison, some error functions A k and A^ are 
displayed simultaneously in Figure [ij One sees clearly that |A fc | > |A fc | on [0, x^], 
while \A k \ < \A k \ on (x k ,oo). 
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Figure 3: Approximation errors A k for k = 0, 1, . . . , 5 (left) and k = 2, 3, . . . , 9 
(right). 



Proof of Theorem [2l Our previous considerations show already that A^ > for 
even k and A k < for odd k on (0, oo). To verify that A k > A k+ 2 for even k and 
A k < A k+ 2 for odd k, we have show that for any integer k > and x > 0, 



(21) 



c* k + (c fe - k)x < x + 



x + 



k + 1 
FT2" 



"fc+2 



+ (cl 



fc+2 



k - 2)x 
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0.0 0.5 1.0 1.5 2.0 

Figure 4: Approximation errors A k and A k for k — 4, 5, . . . , 11. 



But this is equivalent to 

k + 1 " 4 
< =— - x + 



k + 1 



clx + 



- 1 



(22) 



k + 1 - c* 



x + 



C fc+2 + ( C fc+2 ^ 
1 



■ X + 



"fc+2 



according to (12). With 



k + 1- cl 



k + 1 



we may rewrite (|22|) as 
a k x-\ 



1 



(3 k x + 



1 



a k x + 



and 7 fc 



IkX + 1 



C k+2 k 2 



1 



k+2 



IkX+l 



PklkX 2 + /3 k X + 1 

aAlkX 3 + (a k - lk)PkX 2 + (a k + Ik - Pk)% 
PklkX 2 + (3 k x + l 
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But it follows from k + 1/2 < c* k < k + 1 and c* k+2 > k + 5/2 that a k , (3 k , Ik > 0, and 

k + 1 - c* (fc + 2) 2 {k + 1)" 2 4 - k — 2 



& k (* + 2)(* + i)- 

fc + 1 - c* - (k + 2)(/c + l)- 1 ^ + k + 1 
2fc + 2 - (2A; + 3)(A; + 
(fc + 1) 2 
fc + 3/2 



— c 



A' 



= fc + 1/2 + t^— - - c£. 
' 4k + 6 fe 

For fc > 1, 4A; + 6 < 8k + 4 = 8(k + 1/2), so a k - 7 fc > by Lemma [3j while 
«o — 7o ~ 2/3 — 2/7r > 0. Moreover, 

fc + i-4 + (fc + 2)(fc + i)- 1 4-(fe + i) fl 



h = 0. 



fc + 1 



Thus we have verified (21). 

Concerning the comparison of A^ and A k , note that 



^ k + (c k -k)x > ^/c* + (x/2) 2 + x/2 

if, and only if, 



+ (c* k -k- l/2)x > \/c* k + xV4. 
For x > 0, the latter inequality is equivalent to 



2y^(4-fc-l/2) 2^ _ 

* " 1/4 - (c* - fc - 1/2) 2 ~ (cl-k)(k + l-c* k ) Xfc ' 

and we know already that x k is the unique maximizer of A J. On the other hand, 
we also know that the unique maximizer of A J is given by (15), and 



y 'l/4-(cJ-*-l/2)* < 1/2. 
This shows that 

max|A fc (a;)| = |A fc (x fe )| > \A k (x k )\ = |A fc (x fc )| = max|A fc (x)|. 

It remains to show that x k is smaller than one and smaller than the reciprocal 
of y/<%(l - {4:C* k )- 2 ). On the one hand, x = V / V2(4 - ir)/(it - 2) < 1. On the 
other hand, it follows from Lemma|3]that c* k < k + 1/2 + (8c£) _1 , so 

x k < 



1/4 -(8c*)- 2 " v^(l - (4c*)- 2 ) ' 
17 



Note also that the latter bound is strictly decreasing in c* k > 1/2. Since c\ > 3/2, it 
is strictly smaller than 36/35^/2/3 < 0.82 for all k > 1. □ 



6 New Bounds Involving Exponentials 

As a final type of approximation, consider h k as in Lemma [T] with 

g k (x) := x + ^/c* k exp(-5 k x) 
for some 5 k > to be specified later. Here 

g k (x) 2 - xg k (x) — k — g k (x) 

= ^/(f k e-x^(-5 k x){x + ^fd* k ex<p(-5 k x)) - (k + 1) + 5 fc - v /c|exp(— 5 k x) 

k + 1 

~ 5 k + x + a/c£ exp(-<5 fc x) — exp(5 k x) 



=: fk{x)- 
Note that 

/fc(0) = s k - 



k + 1- c* 



( k + 1 

/^.(x) = 1 - 5 k ( y/(% exp(-5 k x) H ^=r exp(5 fc a;) 

V c fc 

/fe(^) = ^(v / 4 ex P(~^ a; ) _ ^-tiexp(5 fc x) 



- 4=(4-^-i) < o. 



k 

Thus /fc is strictly concave on [0, oo) with lim^oo f k (x) = — oo. Setting 

k + l-c 



(23) 4 := 7~^r = V^i 



leads to /jt(0) = and, by Lemma |3j 

c k 

(k + l) 2 



1 + cl- 



> k + 3/2 + (8(k + l))- 1 {kJrl) ' 



6A(k + l) 2 (k + 1/2) + 8(A; + 1) 
18 



k + 1/2 + (8(A; + l))- 1 
> 0. 



Hence, with this choice of 5 k , the function g k satisfies the conditions of Lemma [2j 
Note also that 



[x + ^Jclexp{-5 k x)) - [y^ k + (c* k - k)x) = \fc* k {5 k x + exp(-5 k x) - l) > 0. 
Hence the resulting bounds for 1 — $ are strictly better than the ones in Section [5j 
Theorem 3. Let 

1 



ho(x) := x + v / '2/7rexp(— 5ox), hi(x) := x -\ . _ 

x + v vr/2exp(— 5±x) 

and, for integers k > 2, 

h k (x) := xH , 

x H 



x + 



x + ^Jc* k ey.p(-5 k x) 



where 5 k := ^Jc* k+l — ^Jc\. Then the approximation errors A k := (j)/h k — (1 — $) 
and A k (as in Theorem^ satisfy the following inequalities on (0, oo): 

A > A 2 > A 4 > ■•■ > and A\ < A 3 < A 5 < ■ ■ ■ < 0. 

Moreover, 

| A ft | < |A fe |. 

Figure [5] shows the approximation errors A^ for fc = 0, 2, . . . , 9. Table [l] contains 
some values of the maximum of |A^| and |Afc|, rounded up to four significant digits. 
Although the bounds <fi/h k are much better than (f)/h k or (f)/h k in terms of maximal 
error, note that |Afc(x)| < |Afc(x)| for sufficiently large x > 0. Precisely, |Afe(x)| < 
|Afc(x)| if, and only if, 

exp(d fc a;) > 



and numerical experiments show that this is true for x > 3.2 if k = and x > 3 if 
fc > 1. 

Proof of Theorem [3j We only have to show that A k > A k+2 for even k and 
A k < A k+2 for odd k on (0, 00). This means that for any integer k > and x > 0, 



19 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 

Figure 5: Approximation errors for k = 0,1,..., 5 (left) and k = 2,3, 
(right). 

the difference 

x+ k + 1 Q .( x ) = (k + l)g k+2 (x) , 

X+ k + 2 9k[X) xg k+2 (x) + k + 2 [9k[X) X) 

30 ~ r / x 



9k+2{X) 



(k + l)g k+2 {x) 



c* k exp(-5 k x) 



xg k+2 (x) + k + 2 
~ V c t+i exp(5 k x)g k+2 (x) - xg k+2 (x) - (k + 2) 
(24) = y/c* k+1 exp(5 k x)x - y/c* k+2 exp(-6 k+2 x)x - x' 

+ (k + 2) (exp((5 k - 5 k+2 )x) - 1 



is strictly positive. Note that we utilized (11) twice. Note also that 



S k - 5 k+2 - y/c* k+1 - a/c^ - a/ c* k+3 + a/ c* k+2 

k + 1 , k + 3 



V c t+i + 2 ^ C * k + 2 k + 2 ^°* k+l + V c l+2 

\/ C *k+2 ~ \/ C *k+l _ <Wi 

k+2 ~ k+2 
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7 

k 


max \A k (x)\ 
x>0 1 1 


1 A / \ 1 

max \A k (x) 

x>0 1 1 


1 A { \ 1 

max At i 

X>1 1 1 


1 A / \ 1 

max At x 

x>2 1 1 


1 A ( \\ 

max Afclx) 

x>3 1 1 





2.074 • 10~ 3 


1.571 • 10~ 2 


9.194- lO- 3 


9.374 • lO" 4 


3.550 • 10~ 5 


1 


4.796 • 10~ 4 


3.820 • 10~ 3 


1.606 • 10~ 3 


1.041 • 10~ 4 


2.612 • 10~ 6 


2 


1.723 • 10" 4 


1.622 ■ 10~ 3 


4.687- 10" 4 


1.896 ■ 10" 5 


3.175- 10~ 7 


3 


7.888 • 10~ 5 


8.735 • 10~ 4 


1.764- 10~ 4 


4.591 • lO" 6 


5.226 • 10~ 8 


4 


4.214- HT 5 


5.433 • 10" 4 


7.775 • 10~ 5 


1.342 • 10" 6 


1.059 • 10" 8 


5 


2.499 ■ 10" 5 


3.685 ■ 10" 4 


3.814- 10" 5 


4.480 ■ 10~ 7 


2.497- 10" 9 


6 


1.599- 10- 5 


2.663 • 10~ 4 


2.023 • 10~ 5 


1.655 • 10~ 7 


6.625 • 10~ 10 


7 


1.082 • 10~ 5 


2.010 ■ 10" 4 


1.138 • 10" 5 


6.616 ■ 10" 8 


1.932 ■ 10~ 10 



Table 1: Maximal approximation errors. 



x 



by (12). Hence dividing (24) by a; yields 

v /c|^exp(5 fc x) - ^/c|^exp(-4 +2 x) 

+ (* + 2)(«p(*±!f)-l)/x 
= y^^(exp(5 k x) - 1) - v /c|^ 2 "(exp(-5 fe+2 a;) - l) - x 

> x/Cfc +1 (exp(5 fc x) - l) - y/c* k+2 (exp(-5 k+2 x) - l) - x =: J k (x) 
To verify that J k (x) > for all x > 0, note that J k (0) = and 

4(0) = + ^^2^+2- 1 

= (cl +l -k-l) + (k + 3-c* k+2 )-l 



> 



c k+l c k+2 + 1 

1 1 



l{k + 2) 8(A; + 5/2) 
1 



> 



16(fc + 2)(A; + 5/2) 
by Lemma [3] Finally, 

J 'k\ x ) = V^~i S t ex P^kx) - v /c*^ 2 "^ +2 exp(-5 fc+2 x) 

> V C k+l ~ \/ C l+2 ^k+2 

= ~ k ~ !)4 - (k + 3 - c* k+2 )5 k+ 2 

> 4/2 -W 2 > 0, 
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because c* k+l > k + 3/2 and c* k+2 > k + 5/2. □ 

7 Alternative Representations 

In this section we describe briefly representations of our approximations in terms of 
simple fractions rather than continued fractions. This material is rather standard, 
and we refer to Kouba (2006) for deeper connections with Hermite polynomials. 

Lemma 4. Let ho, hi, h 2 , . . . be constructed as in Lemma[IJ Define polynomials via 
Po(x) := 1, Pi(x) := x, Qo(x) := 0, Qi(x) := 1 and inductively, for k = 2, 3, 4, . . ., 

P k (x) := (k-l)P k _ 2 (x)+xP k ^(x), 
Q k {x) := (k - l)Q k - 2 (x) + xQ k -i{x). 

Then for arbitrary k > 1 and x > 0, 

A;P fc _i(x) + P k (x)g k (x) 



hk(x) 



kQ k -i(x) + Qk(x)g k (x)' 
In particular, if g k (x) = x, then h k (x) = P k +i(x) / Q k+i(x) . 

Table [2] contains a list of the polynomials P k , Q k for < k < 8. Suppose that 
gk(x) = x + Gk{x) for some G k : [0, oo) — » (0, oo), for instance, 

G k (x) := Jet + (x/2) 2 -x/2 = ° h : in Section HI 

1 ^ V k K ' ' ' v/c* + (x/2)2 + x/2 

Gfc(x) := v^ ex p( - (V c fe+i ~~ in Section [6l 



Then 

For instance, 

/i 2 (x) 
/i 3 (x) 
h±{x) 
h 5 (x) 



, Pk+i + PkG k r . 
^fc = 7^ 77777" on[0,oo). 



x 3 + 3x + (a; 2 + l)G 2 (x) 
x 2 + 2 + 2G 2 (x) ' 
x 4 + 6x 2 + 3 + (x 3 + 3x)G 3 (x) 
x 3 + 5x + (x 2 + 2)G 3 (x) ' 
x 5 + 10x 3 + 15x + (x 4 + 6x 2 + 3)G 4 (x) 
x 4 + 9x 2 + 8 + (x 3 + 5x)G 4 (a;) ' 
x 6 + 15x 4 + 45x 2 + 15 + (x 5 + 10x 3 + 15x)G 5 (x) 
x 5 + 14x 3 + 33x + (x 4 + 9x 2 + 8)G 5 (x) 
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k 






P k 


» = 


Ej a J xj 












Qk(x 


) = 


E^- 


X 3 






a 


a\ 


a 2 


a 3 


0,4 0,5 


a 6 






bo 


h 


b 2 






h 


h h 





1 





1 





1 














1 


2 


1 





1 















1 












3 





3 





1 










2 





1 










4 


3 





6 





1 











5 





1 








5 





15 





10 


1 








8 





9 





1 






6 


15 





45 





15 


1 









33 





14 





1 




7 





105 





105 


21 





1 




48 





87 





20 





1 


8 


105 





420 





210 


28 





1 





279 





185 





27 


1 



Table 2: Auxiliary polynomials P k , Q k for < k < 8. 

Proof of Lemma H Note first that h k (x) may be written as 

_ A k (x) + P k {x)g k {x) 
B k (x) + Q k (x)g k {x) 

with certain polynomials A k ,P k ,B k ,Q k . Indeed, 

ho(x) = "tn'^v ^[A (x),P (x),B (x),Q (x)] = [0,1,1,0], 
hi(x) = n~?~* 9l ^\ , so[A 1 (x),P 1 (x),B 1 (x),Q 1 (x)] = [l,x,0,l]. 

For k > 2, let := x + k/g k (x) = (k + xg k (x))/g k (x). Then 

h / x \ _ hk / x \ _ ^fc-i(g)fffc(g) + Pk-i(x)(k + xg k {x)) 
k k ~ l B k _i{x)g k {x) + Q k -i{x){k + xg k (x)) 

kPk-ijx) + (A k _i(x) + xP k -i(x))g k (x) 
kQ k -i(x) + + xQ k -x(x))g k {x) ' 



i.e. 



A fc (x) = kP k -x{x) and P fc (:r) = 4-i(i) + 

B k (x) = kQk-x{x) and = B k -x(x) + xQ k -i(x). 

Since A\{x) = 1 ■ Pq(x) and B\(x) = 1 ■ Qq(x), we may write 

7 / \ kP k ^{x) + P k {x)g k {x) 

hk{x) = —z , s , , , r for fc > 1, 

kQk-^x) + Q k (x)g k (x) 



where 



P k (x) = (k — l)P k - 2 (x) + xP k -i(x) and 

Q k (x) = (k - l)Q k - 2 (x) + xQ k ^i(x) for k = 2, 3, 4,.... □ 
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